clc;clear;

VR = 220e3;
PR = 250e6;
PF = 0.99;

A = 0.9873+0.0071i;
B = 2.85+51i;

%----------------------------------------------%
PR_pahse = PR/3; 
Vphase = VR/sqrt(3);

IR = (PR_pahse)/(Vphase*PF);
current_angle = -acos(PF);

[IR_real,IR_imaj] = pol2cart(current_angle,IR);

VS = A*Vphase + B*(IR_real+i*IR_imaj);

angle = angle(VS) *180/pi